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ABSTRACT 



A theoretical and experimental investigation has been performed 
on the determination of the characteristics of two impinging plane 
turbulent jets for the purpose of understanding the performance of 
proportional fluid amplifiers and momentum-exchange devices. The 
jet deflection angle was determined through the application of the 
free-streamline theory for two normally impinging jets with arbitrary 
throat widths and wall set-backs. The deflection angle and the veloc- 
ity and turbulence distributions across and at various distances along 
the jet were determined through the use of a hot-wire anemometer system. 
The results are presented in terms of normalized parameters suitable for 
comparison with the theoretical predictions. The deflection angles 
predicted theoretically were found to be in good agreement with those 



obtained experimentally. 
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I. INTRODUCTION 



In recent years considerable interest has been shown in the research, 
development, and design of fluid amplifiers with no moving parts. These 
devices are finding many applications in various fields of engineering, 
especially in control systems. The main feature of most of these devices 
is the use of a controlled fluid jet. The power jet or the primary jet 
may or may not attach to one or the other sidewall depending on whether 
the amplifier is a wall-attachment type bistable amplifier or a propor- 
tional amplifier. In the former type, the deflection of the jet from 
one stable to another stable position is accomplished by a combination 
of momentum exchange and pressure build-up in the separation bubble. 

In the proportional or momentum-exchange devices, however, the deflection 
of the jet is accomplished by and proportional to the exchange of momen- 
tum between the primary and control jets. 

Despite the simplicity of their design, the complex combinations of 
fluid dynamic phenomena which take place inside these devices make it very 
difficult to produce a mathematical model able to predict their opera- 
tion. The need for gaining a better understanding of the basic phenomena 
involved in the fluid amplifiers has brought a new interest in the area 
of small flows in closed passages and chambers, and many new works on 
these subjects have appeared in the last few years. 

Some work on the research and development of fluid amplifiers and 
rate sensors has been undertaken by the Department of Mechanical Engineer- 
ing of the Naval Postgraduate School through a research grant from the 
Harry Diamond Laboratories of the U. S. Army Materiel Command. The pre- 
sent work is part of this general research program and deals with only the 



proportional amplifiers . 



A schematic diagram of a typical proportional amplifier is shown 
in Fig. 1. The interaction region of the amplifier is shaped to pre- 
vent the attachment of the supply (power) flow to an output via the 
Coanda effect. With no input (control) flow present, the supply flow 
strikes the splitter and divides equally through the two outputs. If 
pressure is applied to one of the inputs, the resultant control flow 
strikes the supply flow which deflects more of the supply flow to the 
opposite output. The deflection of the supply flow varies proportionally 
with the input flow. The greater the input pressure, the greater is the 
deflection of the supply flow to the opposite output. If the input pres- 
sure differential is reduced, the deflection of the supply flow is re- 
duced proportionately. Since a variation in the control flow conditions 
causes a larger variation at the receivers, the device can be considered 
as an amplifier. 

One can introduce parameters such as differential flow gain, a dif- 
ferential pressure gain, and differential power gain to describe the 
static characteristics of the amplifier. Until recently, these char- 
acteristics were determined purely experimentally without a serious ef- 
fort to theoretically predict the jet deflection, momentum interaction, 
velocity and turbulence distributions, and the division of flow between 
the two outputs and vents. 

In an attempt to lay a theoretical foundation for the prediction of 
the characteristics of proportional amplifiers with one or two control 
flows, the present study was undertaken. Furthermore, in order to partly 
simplify the problem, two two-dimensional, normally-impinging jets were 
considered with different jet widths and side-wall set-backs. Obviously, 
such a configuration (Fig. 2) corresponds to a simple momentum exchange 
proportional amplifier as encountered in several logic devices such as a 
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FIGURE 1. PROPORTIONAL AMPLIFIER 
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1. Auxiliary Jat 

2. Main jat 




PIGURX 2. SCHEMATIC DRAWING OP THE JET ASSEMBLY 




■’half adder" and "transverse impact Modulators" [1]. The understanding 
of the fluid mechanics of such a jet configuration could then be ex- 
tended through the use of similar but perhaps more complicated theoreti- 
cal and experimental means to the study of the complete proportional 
amplifier shown in Fig* 1. 

In connection with this and other current and anticipated investi- 
gations in the laboratory, it was decided to study the impingement of 
two semi -confined jets between two parallel plates, with special regard 
to the influences of the velocity ratios of the two jets and the wall 
set-back. A large scale model was used having dimensions that are large 
compared with the power of resolution of the flow measurement instru- 
ments. Velocity and turbulence measurements were carried out through 
the use of a hot-wire anemometer system for various velocity combina- 
tions and at various sections along the jet. The theoretical analysis, 
based on the free-streamline theory, was used to calculate the jet 
deflection angle as a function of the velocity ratio of the jets assuming 
the fluid to be inviscid. 
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II. THEORETICAL ANALYSIS OF THE JET DEFLECTION 



Many plane irrotational flow patterns which include the formation 
or deflection of free or partly bounded jets can be analyzed satis- 
factorily by the theory of f ree-streamlines by Helmholtz [2] and 
Kirchhoff [3]. The results obtained may be significant not only for 
the corresponding two-dimensional flows of real compressible and in- 
compressible fluids but, in many instances, for their three-dimensional 
counterparts as well. The analysis consists of the definition of suc- 
cessive conformal transformations involving a hodograph, or velocity 
plane, and the application of the Schwartz-Christof f el transformation 
[4]. The fact should be emphasized that from a mathematical point of 
view, the direct calculation of free-streamline flows is restricted to 
two-dimensional irrotational flows of compressible and incompressible 
fluids which are free from gravitational effects. Hence, the forma- 
tion of a jet by two two-dimensional streams impinging obliquely is 
well suited to free-streamline analysis because of the dominance of 
inertia and pressure in the establishment of the resulting flow pat- 
tern and because of the negligible effect of gravity on gas jets. As 
will be discussed later, however, difference in the free-streamline 
velocities of the two jets enhances the shear in the mixing region of 
the two jets. The larger is the velocity difference, the larger the 
energy dissipation becomes. The free-stremaline analysis may be ap- 
plied, however, only when the velocities along both sides of the jet 
are assumed to be identical and equal to an unknown free-streamline 
velocity. Thus the region in which the moment urn -exchange and the 
turbulent mixing take place does not, strictly speaking, follow the 
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assumptions made in the free streamline theory. Hence the validity 
of the analysis depends on the experimental demonstration of the 
uniformity of the velocity profiles at the formation region of the 
resulting iet. 

For the development of the following analysis, the free-streamline 
velocities are assumed to be identical and the validity of the assumption 
is discussed later in connection with the comparison of the theoretical 
results with those obtained experimentally. Mindful of the restrictions 
and limitations of the applicability of this method, the solution of the 
problem under consideration proceeds as follows: 

If 7 and w are used for the complex variable and the complex 
potential, then the variable £ defined as dw/dz is given by 

C = - u + iv - qe - ^ (1) 

in which u and v are the velocity components in the x and y directions, 
n is the magnitude of the velocity vector, and 0 is the angle^ of its 
inclination. 

If the variable Cl is introduced, 




the flow region in the physical plane (Fig. 3a) can be mapped onto a 
corresponding region in the Cl -plane, as indicated in Fig. 3b. Either 
C or 0 must be related to w by means of known methods through still 
another complex variable designated as t. This relationship completes 
the rather complicated functional analysis because, as has already been 
seen, either C or Cl can be expressed in terms of dw/dz. 

Any polygonal boundary in one complex plane can be transformed into 
the entire real axis o f another through application of S 'hwartz-Chris tof fe 1 
theorem, the interior of the polygon being transformed into either the 
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upper or the lower half of the second plane [4]. Thus the outline in 
the -plane is directly transformable into the real axis of the t-plane. 
Therefore, the O-plane shown in Fig. 3b may be transformed to the t- 
plane shown in Fig. 3c. The transforming function can be written as 
follows: 



Upon 



Q = 



Mdt 

v'(t + l) (t-hj 



integration, one has 



(3) 



P = 2M Ln C (v't+1 + v't-h) 



( 4 ) 



in which M and C are found to be 1/2 and 1/ 1 + h by substituting into 
the equation the coordinates of the points B and D. The final form 
of O thus reduces to 



Q = Ln ' / t +l — ±y r td l 



v h+1 

In particular, it will be noticed that at t 
C) and one has 

i£ = Ln — + — = i tan 

Vl+h 



0, Q m i fr 



or 



tan ^ = 7h 



(5) 

(at point 



( 6 ) 



Obviously, the angle of the deflection of the main jet, 0! , is related 
to j9 by the simple relationship 

a + j9 = n/2 (7) 



18 




(a) Physical Plane 




(b) -A - Plane 

p ED C BAP 

rf ~ r o if £ 

(c) t- Plane 



FIGURE 3. DEFINITION OF VARIOUS FLOW PLANES 
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The velocities are related to various parameters in the t-plane 



via the ft transformation. At A, for example, 



ft = 




Ln yfc+1 + s/k-h 

ThTi 



or 



= , = '/hTl 

U j Vk+l + Jk-h 

Likewise, at E 



( 8 ) 



ft = 



Ln 



U m 




Ln 



yr^f + j yF+h 

Jh+l 



or 



jm _ _ y/h+I 

u j ' ~ JT- l + JITh 



( 9 ) 



Before proceeding further. It should be pointed out that the 
parameters (f, h, k) could have been selected at other points or special 
use could have been made of the symmetry of the points B and D. A care- 
ful analysis of the various possibilities has shown that the t-plane 
selected herein will yield the least amount of difficulty in the evalua- 
tion of the integrals to be performed laear. 

The complex potential v can be expressed ln terms of t by the method 
of sources and sinks as 



U4-C U-a U •• 

W » Lnt Ln(t-k) - Ln(t+f) 



( 10 ) 
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From this relationship and from the definition of £ , one obtains 



Z = - 






V t+1 



+ 



a6 



■u c a o 

: " h J It - — 



et) 

t+f 



J dt 



(ID 



The integrals appearing in the above equation may be evaluated for 

various intervals in the t-plane or for various distances between various 

points in the physical plane. The position of the point B relative to 

the point F is given by: 

t=°° 

-t + ia = | z | 

t=h 

which, upon evaluation in a straightforward manner, yields 

| - - — L- - | Ln 4 4 1 + 1 + MTTl Ln 44 + 4+1 
Vh+1 7h+l - 1 v'k+1 - v'h-1 

( 12 ) 

+ ~ -2 T? ^ STl tan 1 - ffT) -f J f+h , 

It should be noted that c/a, 6 , e/a, and T) are related by the equa- 

tion of continuity as 

f = 6 + f n (i3) 

The evaluation of the integrals indicated in Eq. (11) from t - - 1 
to t = - 00 and the separation of the real and imaginary parts similarly 
y ie Ids 

f — ! f it - ff 6/k+l + - rjv'f+h Ln + ^ 1+h 

9 ffv'h+1 " a a v'f+h - v'l+h 

(14) 

- f v'h Ln ^ - 26 v'iTh tan" 1 ” ^ | 

9 v'l+h - v'h V 1+h J 
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It is thus apparent that the velocity ratios and the relative setbacks 
t/a and s/a and the angle jS may be evaluated by assigning various 
values to the parameters f, h, and k in the t-plane. It is further 
obvious that the width of the jet at C may be evaluated without diffi- 
culty by carrying out similar integrations from D to C and from B to C. 
Since the jet width or the contraction coefficient is not needed in the 
present analysis, the results concerning the contraction coefficient will 
not be presented here. 

The expressions for t/a and s/a are of such nature that one cannot 
assign special values to them and expect to find the corresponding 
values of f, g, and k and the angle j3 . The nonlinearity of the equa- 
tions makes it mandatory that one assign arbitrary values to the para- 
meters in the t-plane in a fairly systematic manner and obtain the 
corresponding values of t/a, s/a, and 0 . In this manner, one obtains 
an infinite variety of possible solutions and then selects those which 
conform to the selected geometry. This procedure which is common to all 
free-streamline problems, has been applied to the present problem through 
the use of the IBM 360 computer and the results obtained for both 
t/a=-l 9 s /a * 0, e/a* 1 and for t/a - 0.5 e/a * 1, and s/a * 0 are 
plotted in terms of the jet deflection angle 0L versus the ratio of 
the velocities. These will be discussed later in connection with the 
discussion of the results. The information for other combinations of 
t/a, s/a, and e/a is available in the form of computer print-out and 
they will not be reproduced here for it is practically impossible to 
cover the entire range in the form of a simple plot. Thus only those 
cases which are of practical significance are presented herein and 
compared with those obtained experimentally. The computer program which 
will generate solutions for all ranges of the parameters encountered in 
the analysis is presented in the Computer Output Section. 
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III. EXPERIMENTAL EQUIPMENT 



The experimental apparatus consisted of the jet assembly, supply 
system, hot wire anemometer system, and the velocity calibrator. 

Jet Assembly 

The jet assembly was made of plexiglass for ease of machining and 
flow visualization. As shown in Fig. 2, two ten inch long, one inch 
deep and one-eighth inch wide channels were cut out and sandwiched be- 
tween two one half inch plexiglass plates. The axes of the two jets 
(referred to hereafter as the main ]et and auxiliary jet) were normal to 
each other since in the majority of the momentum-exchange amplifiers 
these two jets are made to intersect each other at an angle of 90 
degrees. In order to prevent premature wal 1 -attachment via Coanda ef- 
fect, the outer edges of the two jets were sharply curved back. The 
three-eighths inch diameter supply lines were connected to the inlets 
of the two channels and a smooth transition was provided between the 
circular and the rectangular sections. In addition, several wire screens 
were placed near the transition region for the purpose of eliminating 
secondary swirling flows which might have developed in the long circular 
supply line. The outer wall of the auxi 1 iary - ]e t channel was made mov- 
able for the purpose of adjusting the relative set-back (see Fig. 2). 

All channel surfaces were polished carefully and the dimensions were 
set with gage blocks to an accuracy in the order of 0.001 M . 

Supply System 

The working fluid (air) was supplied from a 50 Hp compressor. 

Air was then passed through a microfilter for the purpose of removing 
oil, dust, and other impurities and was directed to pressure regulators 
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and Fischer-Porter rotometers. Both the main jet and the auxiliary jet 
were supplied air independently from the main flow line coming out of 
the primary pressure regulator. Each jet had, in addition, a series of 
pressure regulators ranging from ten to 130 psig and an independent flow- 
rator whose capacity was identical with that of the other flowrator. The 
maximum capacity of each flowrator was 14.6 SCFM at 14.7 psia and 70 de- 
grees with a 3/4-27G-10/83 tube and 3/4-GNSVGT-54 float. Finally, in 
order to determine the actual flow rate passing through the flowrators 
at pressures other than 14.7 psia, each flowrator was equipped with a 
pressure gage and the flow rates were obtained by: 

Q * flowrator meter reading x 1007 o full flow xCF(P) 

where 

CF(P) is the pressure correction factor and is provided by Fischer- 
Porter Company in the form of a chart for each flowrator. 

Hot-Wire Anemometer and Calibrator 

The Model 1050 constant temperature hot-wire anemometer of the Thermo- 
Systems, Inc., was used in the measurement of velocities and turbulence 
intensities. The hot-wire probe holder was mounted on a micrometer 
block capable of providing a translation (along a slot along the edges 
of the top and bottom cover plates) and a rotation (in the mid-plane of 
the jets) to the sensor holder. The micrometer enabled the measurement 
of distances along any section in the mid-plane between the top and 
bottom cover plates to an accuracy of 1/160 11 . 

The probe used in the tests had a cold resistance of 6.45 ohms and 
an operating resistance of 9.70 ohms. The probe was calibrated with the 
Model 1125 calibrator of Thermo -Systems , Inc. All calibrations were 
carried out with the high-range calibration orifice of the calibrator 
which had a 0.150" diameter. 
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In order to calculate the velocity of the flow through the orifice, 
the pressure of the calibration chamber was connected to water mano- 
meters of three different ranges to augment the accuracy of the calibra- 
tion. 

A schematic of the experimental layout is shown in Fig. 4. Figures 
5 and 6 are the photographs of the experimental system, the calibrator, 
micrometer assembly, and the plexiglass test model. 
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FIGURE 4. SCHEMATIC DRAWING OF THE EXPERIMENTAL APPARATUS 
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FIGURE 5. GENERAL EXPERIMENTAL SYSTEM, 




Plexiglass jet assembly. 




Model 1125 calibrator. 



FIGURE 6 . JET ASSEMBLY AiJD CALIBRATOR. 
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IV. EXPERIMENTAL PROCEDURE 



The experiments consisted of the setting of proper dimensions to 
the jet channels, the calibration of the hot-wire probe, setting of the 
anemometer, setting of the selected flow rates for each jet, the deter- 
mination of the jet deflection angle a through the use of a flow visuali- 
zation method, and the measurement of velocities and turbulence intensities 
along lines normal to the jet axis in the mid-plane of the jet at various 
distances from the nozzle exit. 

Setting of the Channel Dimensions 
It is apparent that an infinite variety of geometrical configura- 
tions may be tested even if the jet impingement angle is maintained at 
90 degrees. For this reason, attention was given only to those geometries 
which have considerable practical significance as far as the momentum- 
exchange amplifiers are concerned. In the majority of such devices, the 
channel widths are kept identical and the side-wall of the auxiliary jet 
is set back a distance equal to either the nozzle width or one-half of 
the nozzle width. Thus, in the present study, the channel widths were 
kept identical (1/8 ,f ) and two different set-backs were used for the auxi- 
liary jet, i.e., t/a was chosen in the first series of tests equal to 
'0.5 and in the second series of tests equal to unity. All dimensions 
were, as cited before, set with gage blocks to an accuracy of 0.001 M . 
Calibration of the Hot-Wire Probe 

The hot-wire probe was placed 1/16" above the axis of the calibra- 
tion orifice as recommended by the manufacturer [5] and the pressure 
tap of the calibrator was connected to a water manometer. In order to 
increase the accuracy of the pressure measurement, three different 
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manometers were used. These consisted of manometers 0 to 0.2" range 
with 0.01" increments, 0 to 2" range with 0.1 M increments, and a mano- 
meter with 0-36" range with 0.1" increments. The first two manometers 
were of inclined type for further increase of the accuracy in the reading 
of the water elevations. The velocities through the orifice were then 
calculated through the use of the following equation: 



A? is the manometer readings and p the density of air. The 
corresponding voltage output on the hot-wire anemometer was recorded 
for each manometer setting and then a velocity vs. volts "calibration 
curve" was plotted as shown in Fig. 7. 

The calibration curve as obtained above is of limited use because 
of the difficulty it introduces in the evaluation of velocity for each 
anemometer reading corresponding to velocity and turbulence measure- 
ments along various lines, at various distances from the nozzle, and 
for various combinations of the flow rates. It is, therefore, obvious 
that one has to fit one or more polynomials to the calibration data 
and calculate, through the use of the computer, the velocities and 
turbulence intensities corresponding to each anemometer and Root Mean 
Square meter reading. For this purpose, various degrees of polynomials 
were tried for the representation of the entire calibration data and it 
was found that two 4th degree polynomials with a common ordinate and 
tangent at the point of intersection would be most satisfactory. The 
polynomials chosen were of the following type: 



U 



OR 




(15) 



where 




( 16 ) 
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FIGURE 7. CALIBRATION CURVE 



where AE = E - E q and E q is the anemometer output at zero flow velocity 

and U is the temporal mean velocity at a given point- The coefficients 

a of the above polynomial were determined through the use of an IBM 360 
n 

(see Appendix B) . One 4th degree equation with coefficients a was used 

n 

for differential anemometer output equal or below E - E q = 1.460 and 

another 4th degree polynomial with coefficients b^ was used for anemometer 

differential outputs greater than E - E * 1.460. The coefficients a and 

o n 

b are tabulated in Appendix B% 
n 

The calibration curve described above was checked frequently during 
experiments (at least once after the completion of the velocity and 
turbulence measurements for a given flow combination) for the purpose of 
ensuring that neither the probe nor the anemometer system had drifted 
from their original conditions. 

Setting of the Anemometer 

Prior to each experiment, the anemometer was turned on and set to 
warm up for about one-half hour. Then the cold resistance of the hot- 
wire probe was checked to make sure that it had not deviated from its 
previously measured value. Then the operating resistance was set in 
accordance with the regular operational instructions of the anemometer 
[5] and the system was set to 'RUN 1 position. One of the output channels 
of the anemometer was connected to an oscilloscope and the flow was turned 
on and increased gradually until an oscillation, due to feedback in- 
stability in the amplifier circuitries, was encountered. Then the 'TRIM 1 
on the anemometer was adjusted and the oscillation was eliminated. This 
procedure was, of course, carried out for all velocities encountered in 
the experiments by increasing and decreasing the flow rate. 
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The oscilloscope was, of course, observed at all times during the 
actual measurement to make sure that no sys tern-genera ted noise inter- 
fered with the actual turbulence of the flow. 

Setting of the Selected Flow Rates 
The flow rates through each channel were selected on the basis of 
the Reynolds numbers encountered in the application of large momentum- 
exchange devices. The velocity of the main jet was maintained at 173.200 
ft/sec which corresponded to a flowrator reading of 60% of the full capa- 
city and to a Reynolds number of 11300. The flow rates of the auxiliary 
jet ranged from 60% to 107 o at increments of 107 o . These flow rates and 
the corresponding velocities and Reynolds numbers are tabulated below. 

Flow Rate Velocity ft/sec Reynold Number 



60% 


173.200 


11300 


50% 


146.500 


9550 


40% 


113.200 


7380 


30% 


84.096 


5485 


20% 


56.064 


3660 


10% 


28.032 


1832 



During the experiments, a flow rate was selected for the auxiliary 
jet from the above table and the flow through each jet was maintained 
constant (within the limits of the experimental uncertainty) with the 
use of separate pressure regulators described previously. 

The Determination of the Jet Deflection Angle 
The jet deflection angle O' was determined through the use of 
several different procedures. The first two of these procedures con- 
sisted of the use of lamp-black and oil mixture and the use of powdered 
color paint. 
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These gave the jet deflection angle to an accuracy within +2°* 

The third method which yielded the jet deflection angle to a finer 
degree of accuracy was the use of the measured velocity profiles. 

The lamp-black and oil were mixed at more or less arbitrary quanti- 
ties and the inner face of the base plate was carefully painted with a 
brush. Prior to painting, the valves and the pressure regulators were 
adjusted in such a manner that when the main flow was turned on, the 
flow combination would correspond to the desired flow rates in each 
channel. Obviously this latter procedure was easily carried out by 
adjusting the flows prior to the painting and then shutting off the 
main valve rather than the secondary valves. Then when the main valve 
was turned on, the resultant jet partly washed the surface of the paint 
and left a clear print of the jet configuration. It is, of course, 
fully realized that the boundary layer may yield a flow configuration 
on the paint somewhat different from that which would be encountered in 
the mid-plane. Nevertheless, the method yielded the jet deflection 
angles which were fairly satisfactory as far as the preliminary deter- 
mination of the velocity and turbulence traversing planes was concerned. 
The jet deflection angle was measured relative to the main jet by plac- 
ing two graph papers above and below the test section and by lining up 
the paint trace with the lines on the graph papers. 

The second method used in the determination of the jet deflection 
angle consisted of the preparation of a white oil-painted paper and 
the introduction, into the supply lines, of small amounts of powdered, 
colored, dry paint (different color for each channel). The painted paper 
was glued to the inner face of the bottom plate with rubber cement. The 
powdered paint released from the supply line followed the jet configura- 
tion and stuck on the freshly painted paper. Thus, this method quickly 
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yielded the jet configuration and the deflection angle and also gave a 
pretty good idea about the mix-up of the jets in the interaction region. 

In all the experiments cited herein, the mixing region and of course the 
remainder of the jet came out at a color which one would have obtained 
by mixing them thoroughly in a separate container. A sample jet con- 
figuration as obtained in this manner is shown in Fig. 8. 

The Measurement of Velocity and Turbulence Intensities 

The velocities were measured along the lines normal to the iet 
centerline determined through one of the flow visualization techniques 
cited above. The hot-wire probe was moved at suitable increments along 
lines normal to the jet centerline and both the velocity and Root Mean 
Square (RMS) meters were read and recorded. The velocities and, of 
course, turbulence intensities were measured at six different stations 
for six different flow combinations. These stations were located at 0.15", 
0.3 M , 0.8", 1.5", 2.5", and 4" from the point where the inner walls of 
the two iets intersected. As will be mentioned later, only 34 velocity 
and 34 turbulence intensity traces are presented instead of 36 of each. 
(Because the flow rates of 1/6 and 2/6 made it impossible to maneuver the 
hot-wire probe at the first station (0.15 M from the nozzle) due to the 
fact that the traversing line fell within the auxiliary iet channel.) 

Figure 9 is a schematic drawing of the velocity traversing line rela- 
tive to the jet channels. 

The turbulence intensity 



/ - 2 /_ 
V u / u 



was calculated by 



vu 2 /- = M 

/u u dl 



(17) 



2 



where Ve is the output of the RMS meter and dU/dE is the deriva- 
tive of Eq. 16 or, in other words, the complement of the slope of the 
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FIGURE 9. SCHEMATIC OF VELOCITT TRAVERSING LINE 
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calibration curve. The derivation of the above equation is shown in 
Appendix C. The turbulence intensity profiles were plotted like the 
velocity profiles by means of an IBM 360 computer and will be discussed 
subsequently. 
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V. DISCUSSION OF RESULTS AND CONCLUSIONS 



The purpose of the present investigation was the understanding of 
the characteristics of the impingement of two plane jets. These character- 
istics essentially were the deflection angle of the resulting jet, and the 
velocity and turbulence distributions. 

The deflection angles obtained experimentally are shown in Fig. 10 
for a 1/16 M set-back (t/a = -0.5). The similar data for a 1/8 M set-back 
(t/a = -1.0) are shown in Fig. 11. Also shown on both figures are the 
results obtained theoretically through the use of Eqs . (12) and (14). 

It is apparent from the examination of Figs. 10 and 11 that the 
agreement between the theoretical and experimental results is quite 
satisfactory, at least for the flow rates and set-backs used. 

The velocity distributions for various angles of deflection at 
various distances along the nozzle exit are shown in Figs. 12-29 in a 

normalized form, i.e., U/U versus x/a. The following facts become ap- 

m 

parent upon a brief examination of these figures: 

(a) The velocity profiles become flatter and perfectly sym- 
metrical as the distance increases from the nozzle exit to the section 
where the velocity is measured. This, of course, is an observation 
which is true for all jet flows. 

(b) The mixing of the two jets takes place more quickly for 

the smallest jet velocity ratios such as U /U = 1/6. As the velocity 

a m 

ratio increases, the velocity profile exhibits two maximums with a 
trough in the middle. This points out that the mixing of high velocitv 
jets requires larger distances along the jet but certainly a distance 
less that approximately d/a = 6. The significance of this observation 



39 




o o o o o o 

in Tt m cn pH o 






CN 



00 

o 



KO 

o 



o 



CN 

o 



o 

o 




« 



40 



FIGURE 10. JET DEFLECTION ANGT.F 
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FIGURE 11. JET DEFLECTION ANGLE 



is in relation with the theoretical analysis. As has been cited in 
the development of the theoretical analysis, the f ree -s treaml ine velo- 
cities on the two sides of the inviscid jet had to be assumed to be 
identical. (The velocity profiles for all velocity ratios, near the 
nozzle where the free -streamline analysis is supposed to hold true, show 
an almost uniform distribution with the exception of two peaks and a 
trough whose amplitude is about 1/8 of the maximum of the normalized 
velocity. The velocity profiles for each angle are combined in Figs. 
30-35 for the purpose of giving a quick idea about the evolution of the 
jet from the nozzle downstream. 

Figures 30-41 are plotted using the parameters U/U^ vs. x/d. The 
purpose of these plots was to explore the similarity in the velocity 
profiles in the initial region of jet development. Previous studies [6] 
on the evaluation of single two-dimensional and axi -symmetric jets have 
shown that plots prepared in the manner cited above resulted in a single 
normalized curve. This was, to be sure, no mixing occurs with the ex- 
ception of that with the surroundings. The plots presented herein show 
that such a similarity exists only for distance of d/a > 2.4. This is 

evidenced by the fact that the two velocity profiles at distances of 
d/a * 1.2 and d/a 28 2.4 do not fall on the rest of the velocity profiles 
exhibiting similarity. This may be interpreted as a consequence of the 
mixing of two jets of unequal velocity. 

The turbulence intensities were plotted as a function of x/a for 
various angles of delfection and for various stations along the jet and 
are shown in Figs. 42-59. The main characteristic of all the turbulence 
distributions is that there are two high regions of turbulence and a 
region of relatively low turbulence in the central core of the jet. The 
result is not unexpected in view of the fact that aach Jet has the 
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highest shear and mixing with the surrounding ambient flow at its two 
sides. It is also apparent from these figures that the intensity of 
turbulence increases with increasing distances along the jet. This is 
not so much because of the absolute increase of the RMS value of the 
local values of velocity fluctuations, but rather because of the rela- 
tive decrease of the magnitude of velocity along the jet. The turbulance 
intensity profiles for each angle are combined in Figs. 60-65. 
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PIGURE 12 NORMALIZED VELOCITY PROPILES 
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FIGURE 13 NORMALIZED VELOCITY PROFILES 
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PICOU lif NGBULIZED VILOCITT PKOPILES 
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FIGURE 15 NORMALIZED VELOCITY PROTILES 
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FIGURE 17 NORMALIZED VELOCITY PROFILES 
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FIGURE 19 NORMALIZED VELOCITY PROFILES 
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FIGURE 21 NORMALIZED VELOCITY PROFILES 
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FIGURE po NORMALIZED VELOCITY PROTILES 
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FIGURE 2k NORMALIZED VELOCITY PROFILES 
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FIGURE 25 NORMALIZED VELOCITY PROFILES 
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PIGUIX 26 NOKMtLIZED TELOCITY PMPILKS 
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FIGURE 27 NORMALIZED VELOCITY PROFILES 
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FICU1K 28 H0KML1ZED VKLOCITT FKOriLSS 
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FIGURE 29 NORMALIZED VELOCITY PROFILES 
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PIGUKX on WOtMLLIZZD VELOCITY PKOPILXS 
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FIGURE 31 NORMALIZED VELOCITY PROFILES 
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LIZES VELOCITY PSOPILZS 




SO 



90 

* 



CM 

co 



o 

o 



sO 

• 

H 

I 



CM 

CO 



00 

t 



«» 

• 

vO 

I 



m 

K 



to 



65 



FIGURE 33 NORMALIZED VELOCITY PROFILES 
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FIGURE U NQKMALI2ED VELOCITY PROFILES 
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FIGURE 35 NORMALIZED VELOCITY PROFILES 
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FIGURE 36 NORMALIZED VELOCITY PROFILES 
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FIGURE n 7 NORMALIZED VELOCITY PROFILES 
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FIGURE -jQ NORMALIZED VELOCITY PROFILES 
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FIGURE 39 NORMALIZED VELOCITY PROFILES 
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FIGURE l^o NORMALIZED VELOCITY PROFILES 
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FIGURE jj/i NORMALIZED VELOCITY PROFILES 
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FIGURE h2 TURBULENCE INTENSITY PROFILES 
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FIGURE 43 TURBUIANCE INTENSITY PROFILES 
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FIGURE U TURBULENCE INTENSITT PROFILES 
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FIGURE | r' TURBU1ANCE INTENSITY PROFILES 
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FIGURE i/ TURBOLAHCE INTENSITY PROFILES 





79 



FIGURE lj.7 TURBULENCE INTENSITY PROFILES 
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FIGURE 4 9 TURBULENCE INTENSITY PROFILES 
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FIGURE c:o turbulence intensity profiles 
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FIGURE 51 TURBULANCE INTENSITY PROFILES 
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FIGURE co TUtBOLANCE WINS ITT PROFILES 
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FIGURE 33 TURBULENCE INTENSITY PROFILES 
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FIGURE 5^. TURBULANCE INTENSITY PROFILES 
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FIGURE TURBULANCE INTENSITY PROPILES 



Ua/Ura ( 56.0/173.2) 



>6 



O 

O 

It 



fa 

c 



o 



s 

-p 



-d* o 
cvj o o 

P^vO ^ 

• vO T* O 
O -r- 

tl II II II 



O 

•P -P H 
O O 50 

* r ~3 £j 




<0 © 



vO 

o 



\r> 

O 







88 



FICUU tji TURB01ANCE HfTKMSITY PR Of I LIS 
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FIGURE 57 TURBULENCE INTENSITY PROFILES 



vP 




(D 

3 



|"tg 



90 



FIGURE 58 TURBULANCE INTENSITY PROFILES 
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FIGURE ^9 TURBUIANCE INTENSITY PROFILES 
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FIGURE £q TURBULANCE INTENSITT PROFILES 



o 




sO 

o 



• • 
o o 



cn 

o 



<NI 

O 



o 

o 



00 

<± 



CM 

m 



O 

O 



sO 

I 



<N 

• 

cn 

i 



00 

t 



<* 

• 

sO 

I 



« 

X 



ID 




93 



FIGURE TURBUIANCE INTENSITY PROFILES 
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FIGURE 42 TUUOLAMCZ ivtensity profiles 
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FIGURE TURBULANCE INTENSITY PROFILES 
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FIOJKE Q. TURWJLAHCE I *T IMS ITT PROTILES 
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FIGURE 65 TURBULANCE INTENSITY PROFILES 



V I . RECOMMENDATIONS FOR FURTHER WORK 



The following are recommendations for further studies: 

a. Jet deflection, velocity and turbulance distributions be 
studied for other set-back and channel widths. 

b. The study of the above parameters be extended to lower 
Reynolds numbers where one or both of the jets are laminar. 

c. A theoretical analysis be carried out for unidentical 
free-stream velocities, and finally, 

d. The aspect ratio (the ratio of the width of the channel 
to the depth of the channel) be changed from 2 to 8. 
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APPENDIX A 



Evaluation of the Integrals 

The integrals indicated by Eq. 11 may be evaluated as follows: 



I 



1 



at = c { 2/£Tl - Ln aft 71 + } 

yt+i - yr 



I 



2 



a 6 



’ ^1 
J t-k 



dt 



b 6 { 2JW. - ,/kJI Ln -ft* 1 * 

Jt+1 - Jk+i 



I 



3 




•J t +1 
t+f 



dt 




2.\/t+i + 2\/f-T 




} 



I A = c 



t-h 

t 



dt = c | 2j t-h + 2 \/h tan 



-1 



t-h 



I 



5 



a 6 



/t^h 

t-k 



dt 



-a 6 {2/t=h - ^/k^h Ln + } 

v't^h - v/k^h 



I* = - er? 



dt = -er) { 2 ,/t^h + 2 v^f+h tan 1 J - } 
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The above results are evaluated between proper limits to yield the 
normalized distances s/a and t/a given by Eqs. (12) and (14). Note 
that regardless of the integration limits the sum of the first terms 
in the integrated expressions is zero because of continuity (see Eq. 
(13)). 
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APPENDIX B 



Determination of the Coefficients of the Polynomial 
One of the computer programs for the determination of the coeffici- 
ents of the polynomials describing the two parts of the calibration curve 
is shown on page 

The polynomials for either part of the calibration curve were as- 
sumed to be of the following type. 

U(E) = a 1 + a 2 AE + a 3 (AE) 2 + a 4 (Ae) 3 + a 5 (AE) 4 

The coefficients a and b for the two polynomials are tabulated below: 
n n 





n 




n 


a l 


= 0.764 


b l 


= 3728.036 


a 2 


= 17.491 


b 2 


—9160.219 


a 3 


= 20.558 


b 3 


= 8384.484 


a 4 


=-20.184 


b 4 


—3353.178 


a 5 


= 24.969 


b 5 


= 512.192 



The second program used in the determination of the coefficients is 
not presented herein for it was identical in every aspect with the pro- 
gram used for the first part of the calibration curve with the excep- 
tion of course, of the test data and the range of iteration. The 
resulting calibration curve and the test data are shown in Fig. 7. 
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APPENDIX C 



Analysis of the Turbulence Measurement 



Evaluation of the intensity of the turbulence is based on the well- 

2 2 Vh 

n King’s Law given as E -E q 31 BU where 

E is the anemometer output at zero velocity, E is the output 



for the velocity U, and B is a constant. The exponent n is ordinarily 
assumed to be 2. In most of the hot-wire work covering a large range 
of velocity, the exponent may deviate considerably from 2 and it is 
for this reason that in precision work one must use a calibration curve 
The evaluation of such a calibration curve was described in the pre- 
vious appendix. In the following the evaluation of the corresponding 
turbulence intensities will be described. 

For a turbulent flow the instantaneous value of the anemometer out 



put is given by E = _ 



Where E is the voltage output corresponding to the mean velocity at 



o 




a given point and 




is the Root Mean Square value of the 



voltage output corresponding to 




(RMS value of the fluctuating 



component of velocity). 



Writing 



E = f(U) 



i 



dE = f ' (U) dU 



and 




or 



f (U) 
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and noting that 




Thus the intensity of turbulence reduces to 




U Uf'(U) U dE 



where 

f 1 (U) 



dE 

dU 



The turbulence intensity 




was evaluated once again through 



the use of computer 
and the values of 
respective parts of 



by feeding in the measured values of 
dU/dE* aS °b ta i nec * from the two 
the calibration curve. 




polynomials for 
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